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ABSTRACT 

Context. The forthcoming arrival of the NASA's Dawn spacecraft to Ceres and Vesta means that these unexplored worlds in the 
asteroid Main Belt are targets of several studies. 

Aims. We study the impactor flux and cratering on Ceres and Vesta caused by the collisional and dynamical evolution of the asteroid 
Main Belt. 

Methods. We develop a statistical code based on a well-tested model for the simultaneous evolution of the Main Belt and NEA size 
distributions. This code includes catastrophic collisions and noncollisional removal processes such as the Yarkovsky effect and the 
orbital resonances. It is worth noting that the model assumes that the dynamical depletion of the early Main Belt was very strong, and 
owing to that, most Main Belt comminution occurred when its dynamical structure was similar to the present one. 
Results. Our results indicate that the number of Z) > 1 km Main Belt asteroids striking Ceres and Vesta over the Solar System history 
are approximately 4 600 and 1 100, respectively. Moreover, the largest Main Belt asteroids expected to have impacted Ceres and Vesta 
had a diameter of 7 1 .7 for Ceres and 21.1 km for Vesta. As for the cratering, our simulations show that the surfaces of Ceres and Vesta 
present a wide variety of craters with different sizes. In fact, the number of Z) > 0. 1 km craters on Ceres is ~ 3.4 x 10^ and 6.2 x 10' 
on Vesta. Moreover, the number of craters with D > 100 km are 47 on Ceres and 8 on Vesta. On the other hand, our study indicates 
that the D = 460 km crater observed on Vesta had to be formed by the impact of a £) ~ 66.2 km projectile, which has a probability of 
occurr ~ 30 % over the Solar System history. 

Conclusions. If significant discrepancies between our results about the cratering on Ceres and Vesta and data obtained from the Dawn 
Mission were found, they should be linked to a higher degree of collisional evolution during the early Main Belt and/or the existence 
of the late heavy bombardment. An increase in the collisional activity in the early phase may be provided for an initial configuration 
of the giant planets consistent with, for example, the Nice model. From this, the Dawn Mission would be able to give us clues about 
the initial configuration of the early Solar System and its subsequent dynamical evolution. 
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1. Introduction 

Ceres and Vesta are the largest and most massive members of 
, a vast population of small bodies located between Mars and 

■ Jupiter commonly called the asteroid Main Belt. Several nu- 
; merical models (e.g. Bottke et al. I2005al l2005bl O'Brien & 
, Greenberg 2005; de Ella & Brunini l2007l l indicate that the size 

distribution of the Main Belt asteroids is determined primarily 
by collisional processes. In fact, these studies suggest that most 
, of the largest objects (D > 120 km) have never been disrupted, 

■ while many smaller asteroids are byproducts of fragmentation 
events among the largest bodies. While Ceres and Vesta have 
not been targets of catastrophic collisions, they have been ex- 
posed to cratering impacts over the age of the Solar System. In 
fact, the existence of the Vesta family (Binzel & Xu 1 19931 1 is 
clear evidence that this object has undergone large cratering im- 
pacts over time. Cratering is one of the most important processes 
that determine the morphology of the surface of a Solar System 
object. The understanding and quantification of the impactor 
source population onto an object and the observation of the ob- 
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ject surface help for understanding the dynamical and physical 
history of both the impactor population and the target. 

Launched in September 2007, NASA's Dawn Mission was 
captured in orbit by Vesta on July 15, 201 1, and it should reach 
the vicinities of Ceres in February 2015. The theoretical predic- 
tions of producing of craters may be compared with observations 
of Ceres and Vesta. This will help, on the one hand, for identify- 
ing the source of craters and, on the other hand, accounting for 
the geological processes that have acted on the surfaces of those 
bodies. Then, it is very important to study all the possible sources 
of crater production on Ceres and Vesta in order to estimate the 
total crater production and to contrast them with observations. 

In this paper we evaluate the impactor flux and cratering on 
Vesta and Ceres due to the collisional and dynamical evolution 
of the asteroid Main Belt. To do this, we constructed a statisti- 
cal code based on the collisional model developed by Bottke et 
al. (I2005al) with some dynamical considerations from Bottke et 
al. (I2005bl) . A comparison between our study and data obtained 
from the Dawn Mission may be relevant for the structure and 
evolution of the early Solar System. 
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2. The full model 



Bottke et al. (I2005al l developed a collisional model capable of 
tracking the evolution of the asteroid Main Belt over the Solar 
System history. At each timestep, this statistical algorithm calcu- 
lates the total number of catastrophic collisions between objects 
residing in different size bins, using parameters such as the mean 
impact velocity (V), the intrinsic collision probability (Pi), and 
the impact energy required for dispersal Qo- From this, the code 
computes how many objects are catastrophically fragmented and 
removed from each size bin, as well as the number of fragments 
resulting from those collisions, which are distributed in the dif- 
ferent bins according to their sizes. This model assumes that all 
breakups occur close to the catastrophic disruption threshold and 
neglects cratering events, which produce much less ejecta over 
time than catastrophic disruption events, and highly-energetic 
catastrophic disruption events, which are relatively uncommon. 

The algorithm developed by Bottke et al. (2005a) includes 
the effects of an intense period of collisional evolution in the 
early massive Main Belt. Using numerical simulations, the au- 
thors found that the net collisional activity in the Main Belt over 
its lifetime is the equivalent of ~ 7.5 - 9.5 Gyr of collisional 
activity in the current Main Belt. This "pseudo-time approxi- 
mation" means that the Main Belt population required ~ 1.5-2 
times the degree of comminution than it would have experienced 
if it were not initially much more massive. It is worth noting that 
this model assumes that the collisional history of the Main Belt 
has been dominated by the same intrinsic collision probabili- 
ties and impact velocities as are found in the current Main Belt. 
Based on Petit et al. (i200 II 120021) . Bottke et al. (2005a) consider 
that the dynamical removal phase of the Main Belt was short 
and owing to that the high velocity impacts did not play an im- 
portant role in the collisional history of the Main Belt. From this, 
the collisional model from Bottke et al. ( I2005al ) uses the current 
values of {Pi) and (V) throughout the whole simulation. 

This model is able to fit a wide range of observational con- 
straints such as the Main Belt size distribution, the number of 
large asteroid families produced by the disruption of D > 100 
km parent bodies over the past 3 to 4 Gyr, the existence of a D = 
460 km crater on the intact basaltic crust of Vesta, and the rela- 
tively constant crater production rate of the Earth and Moon over 
the last 3 Gyr. 

Later, Bottke et al. ('2005b') performed a study aimed at link- 
ing the collisional history of the asteroid Main Belt to its dy- 
namical excitation and depletion. This work combines the colli- 
sional evolution code created by Bottke et al. (2005a) with dy- 
namical results from Petit et al. ( 120011 ). as well as the removal of 
bodies from the Main Belt due to the action of resonances and 
the Yarkovsky effect, which enter the near-Earth asteroid (NEA) 
population. This collisional and dynamical evolution code also 
satisfies the above constraints and successfully reproduces the 
observed NEA size distribution. It is worth noting that Bottke et 
al. ( I2005bl ) validate the pseudo-time approximation proposed by 
Bottke et al. (I2005al ). 

2.1. Simulation parameters 

Here, we construct a colUsional model based on the one de- 
scribed in Bottke et al. ( 12005 al ) with some dynamical consider- 
ations taken from Bottke et al. (p005b ). From this, we track the 
simultaneous evolution of both the NEA and Main Belt popula- 
tions by simulating the effects of an intense collisional evolution 
in the early massive Main Belt from the pseudo-time approxima- 
tion proposed by those authors. 
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Fig. 1. Initial Main Belt size distribution used in our simula- 
tions, the observed Main Belt size distribution from Bottke et 
al. ( 2005a> . which is based on Jedicke et al. (2002) with a few 
changes (see Sect. 2.2), and the observed NEA size distribution 
from Stuart & Binzel (|2004| ). 



The initial population used here follows the idea proposed 
by Bottke et al. (I2005all2005bl) . In fact, our starting size distri- 
bution for D > 200 km uses the number of observed Main Belt 
asteroids, with a few objects added in, to account for the Eos and 
Themis parent bodies. For ~ 120 < D < 200 km, the population 
follows an incremental power-law index of -4.5, while for D < 
120 km an incremental power-law index of -1.2 is assigned. On 
the other hand, the NEA population always starts with no bod- 
ies. Figure[T]shows the initial Main Belt size distribution used in 
our simulations, the observed Main Belt size distribution from 
Bottke et al. ( I2005al ). which is based on Jedicke et al. (20021) 
with a few changes, and the observed NEA size distribution from 
Stuart & Binzel (2004). Jedicke et al. (2002 ) analyzed the abso- 
lute magnitude H distribution of the Main Belt in the range 5.0 
< H < 18.5. To transform the Jedick e et al. (120021) H distribution 
into a size distribution, Bottke et al. (12005 al ) assume a visual ge- 
ometric albedo = 0.092. Moreover, they include the observed 
asteroids for D > 300 km using the IRAS/color-albedo-derived 
diameters cited in Farinella & Davis ( 19921 ). 

For collisions between Main Belt asteroids, the value 
adopted in the collisional algorithm for the intrinsic collision 
probability is (Pi) = 2.86x 10 km ~^ yr and the mean 
impact velocity is (V) - 5.3 km s"' (Bottke et al. 11 9941 ). To 
compute impacts on Ceres and Vesta, we use the particular val- 
ues of (Pi) derived by O'Brien et al. (1201 11 1 for each of these 
bodies. The values adopted for (Pi) are 3.70 x 10"'^ km"^ yr"' 
for Ceres and 2.97 x 10""* km yr "' for Vesta. As for the 
disruption law, we use the impact energy required for dispersal 
Qo derived by Benz & Asphaug (1999) for basalt at 3 km s 
While the mean impact velocity (V) for collisions between Main 
Belt asteroids is 5.3 km s"', we adopt the Qo law for 3 km s"' 
since Bottke et al. (2005a) show that it produces very good fits 
to the observational constraints. 

The differential fragment size distribution (FSD) produced 
by each catastrophic disruption event is represented by 



dN = BD-i'dD, 



(1) 



where D is the diameter, dA^ is the number of fragments in the 
size range (D,D + dD), Bis a constant, and p is the power-law in- 
dex. The FSDs used here are based on those observed in asteroid 
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families Uke Themis or Flora (Bottke et al. I2005al) . According 
to Tanga et al. ('1999'), the Themis family was produced by the 
super-catastrophic breakup of a D = 369 km body, while the 
Flora family was created by the barely-catastrophic breakup of a 
D - 164 km body. Bottke et al. (2005a) show thatFSDs resulting 
from super-catastrophic disruption events are represented well 
by a single p value, although this is not true for FSDs produced 
by barely-catastrophic breakups. From this, Bottke et al. (2005a]l 
developed two different FSDs to describe the possible outcomes 
of a catastrophic collision. On the one hand, for D > 150 km dis- 
ruption events, the diameter of the largest remnant is assumed to 
be 50 % the diameter of the parent body. Moreover, the power- 
law index p of the differential FSD between the largest remnant 
and 1/60 the diameter of the parent body is -3.5, while the p 
value is -1.5 for smaller fragments. On the other hand, for dis- 
ruption events among D < 150 km bodies, the largest remnant is 
80 % the diameter of the parent body. Moreover, the power-law 
index p of the differential FSD between the largest remnant and 
1/3 the diameter of the parent body is -2.3. In this case, the FSD 
has a second break at 1/40 the diameter of the parent body. The 
p value between 1/3 and 1/40 the diameter of the parent body is 
-4, while p is assumed to be -2 for smaller fragments. 

On the other hand, since the NEA population is several or- 
ders of magnitude smaller than the Main Belt population, we 
neglect collisions between NEAs and between NEAs and Main 
Belt asteroids. In fact, we consider that the NEA population is 
sustained by the input of material from the Main Belt via or- 
bital resonances and the Yarkovsky effect. Here, we use the non- 
collisional removal rate of asteroids from the Main Belt pro- 
posed by Bottke et al. (2005b). Thus, bodies removed from the 
Main Belt by noncollisional processes are placed in the NEA 
population, which is assumed to decay with a mean dynamical 
lifetime of 5 Myr. 

2.2. Runs 

We perform our numerical simulations using a pseudo-time of 
9.5 Gyr In addition, since the range of 7.5 to 9.5 Gyr is what is 
found by Bottke et al. (2005a) to give the best fits, we decide to 
carry out numerical simulations for a pseudo-time of 7.5 Gyr in 
order to test the sensitivity of our results to this parameter. This 
point is referred in Sect. 3. 

On the other hand, to simulate the disruption events in the as- 
teroid Main Belt in a more realistic manner, our algorithm works 
stochastically treating the breakups as Poisson random events 
(Press et al. 1989 ). A stochastic code produces different results 
using different seeds for the random number generator. Thus, we 
develop a large number of runs using different random seeds and 
then interpret the results statistically. 

To obtain a quantitative measure of how good a run repro- 
duces observational data, we follow the procedure described by 
Bottke et al. (2005a). First, the metric used to determine the 
goodness of fit between the observed Main Belt size distribution 
(Nmb) and the model results (A^model) is given by 



Table 1. Main Belt model parameters (Bottke et al. l2005al l. 
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I Nmodel(D) - Nmb(D)\ 



0.2Nmb(D) 
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where the summation extends over 1 < D < 1 000 km size bins. 
Second, a standard x^am ^^^^ used to compare the number of 
breakups produced in each D > 100 km size bin over the last 
3.5 Gyr to the observed number of asteroid families in the 100 
< D < 400 km size bins. Table[T]summarizes the observed Main 
Belt model parameters necessary to develop these tests (Bottke 
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H is the absolute magnitude, D the central diameter (in kilometers) of 
the bin for a visual geometric albedo = 0.092, dN the incremental 
number of asteroids in each bin based on Jedicke et al. ( 2002 ) with a 
few modifications (see Sect. 2.2), and rfA'pam the number of observed 
asteroid families in each bin. 



et al. l2005at . Following Bottke et al. ( I2005bl l. we assume that a 
run produces a positive match if (Asfd ^ ^'^'^^fam ^ 
probability > 5 %). 

Given the stochastic nature of our code, we carried out 
10000 numerical simulations using a pseudo-time of 9.5 Gyr, 
533 of which show positive matches. From such simulations, 
we calculated the flux of Main Belt asteroids of different sizes 
onto Ceres and Vesta. Besides using a suitable expression we 
also computed the crater diameters produced by Main Belt as- 
teroids on such bodies. Below, we discuss the algorithms that 
relate crater diameter to impactor diameter. 



2.3. Cratering laws on Ceres and Vesta 

Cratering processes have been extensively study through impact 
and explosion experiments (e.g. Schmidt & Housen 119871) . In 
a recent work interpreting the observations of the Deep Impact 
event, Holsapple & Housen (2007 ) have updated the crater scal- 
ing laws for different materials. The impact cratering scaling 
laws depend in general on two exponents yU and v and a constant 
Ki that characterize the different materials. The diameter Ds of 
a crater produced by an impactor of diameter d can be obtained 



4 



G. C. de Elia & R. P. Di Sisto: Impactor flux and cratering on Ceres and Vesta 



Table 2. Values of the mass M, the density p (Baer & Chesley 
|2008] |. the diameter D (Thomas et al. ll997bll2U05T l. and the sur- 
face gravity g used for Ceres and Vesta. 
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from Holsapple & Housen (12007! ) from the equation: 
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(3) 

where p, is the target density, g its superficial gravity, Y its 
strength, p, the density of the impactor, and v,- the impactor ve- 
locity. 

To calculate the craters on Ceres and Vesta, we have to 
take their composition and superficial characteristics into ac- 
count. Table |2] shows the values of several physical parameters 
for Ceres and Vesta. We have no direct information on the sur- 
faces of those asteroids. For Ceres, observations from their re- 
flectance spectra plus theoretical analysis are consistent with a 
rocky surface (Zolotov |2009l l. The low density of Ceres implies 
that it contains low-density compounds. The internal structure 
and the amounts of water in Ceres have already been modeled 
and studied (McCord & Sotin 120051 Zolotov [20091 1, but there is 
no unique, generally accepted model for Ceres interior Those 
models may really be tested with the observations of the Dawn 
Mission in 2015. Vesta has a basaltic surface and is considered 
a differentiated body (McCord et al. 1 19701 1. Then, Ceres and 
Vesta can be assimilated to the cratering law that corresponds 
to wet soils and rock. For this material fi = 0.55, v - 0.4, and 
Ki - 0.93. Then, replacing these values in Eq. Q we can obtain 
the final diameter of the crater for a given impactor diameter. 

Equation (O gives a general law for impacts, and it is a 
convenient empirical smoothing function to span the transition 
between the gravity regime and the strength regime (Holsapple 
I1993I I. We use this general form, since in the low gravity of aster- 
oids the strength regime can be important for the smaller craters. 
The strength of the target F is a value that have many measures 
for a geological material. Housen & Holsapple (120031 ) show that 
the crater size law depends on the average of the tensile and com- 
pressive strength. Then, we adopt a value of the strength Y cor- 
responding to soft rocks of 3 MPa as in Holsapple & Housen 
(l2U07l l. 

Obtained from Eq. ([3]), D, is the diameter of a simple crater, 
but for a large crater, gravitational forces modify the crater form 
leading to complex craters. The transition from simple to com- 
plex craters is an observable quantity unique for each world. 
Schenk et al. (^2004) summarize results for the single-to-complex 
transition diameters for silicate planets and icy satellites. From 
this, an inverse correlation between transition diameter and sur- 
face gravity is evident, with icy satellites offset toward lower 
transition diameters. Since the surface gravities of Ceres and 
Vesta are comparable (Table |2]), extrapolation of the values as- 
sociated with the silicate planets leads to a transition diameter 
D, for those objects of ~ 50 km. When we use the diameter of 
the complex crater from McKinnon et al. ( 120031) . then the final 




E 



E 
o 



1e+08 
1e+07 
1e+06 
100000 
10000 
1000 
100 
10 
1 



Jedicke et al. (2002) 
Final Main Belt Population 
Stuart & Binzel (2004) 
Final NEA Population 



0.1 



10 

Diameter (km) 



100 



1000 



Fig. 2. Our estimates of the Main Belt (solid line) and NEA 
(dashed line) size distributions for a pseudo-time of 9.5 Gyr 
Moreover, the observed Main Belt size distribution from Bottke 
et al. ( I2005al l. which is based on Jedicke et al. ( 120021 ) with a few 
changes, and the observed NEA size distribution from Stuart & 
Binzel (i2004j are also shown. 



crater diameter is given by 
D = D, 

D = LllDADJDtf" 



for D, < D,, 
for Ds> D,. 



(4) 



3. Results 

We present here our main results for the impactor flux and crater- 
ing onto Ceres and Vesta caused by the collisional and dynamical 
evolution of the asteroid Main Belt. 

The statistical analysis developed in this section is based on 
the 533 runs that produce positive matches with the observa- 
tional data. In fact, Figure|2]shows our estimates of the Main Belt 
and NEA cumulative size distributions for a pseudo-time of 9.5 
Gyr, obtained from the 533 runs mentioned above. The results 
obtained from the simulations using a pseudo-time of 7.5 Gyr 
do not show any significant changes in the Main Belt and NEA 
size distributions and the formation of the large asteroid fami- 
lies. As for the impactor flux on Ceres and Vesta, the number of 
impacts on each of these bodies estimated using a pseudo-time 
of 7.5 Gyr is 70 to 80 % what is obtained with a pseudo-time of 
9.5 Gyr 



3.1. Impactor flux onto Ceres and Vesta 

Figure |3] shows the total number of Main Belt asteroid impacts 
onto Ceres and Vesta as a function of impactor diameter obtained 
using a pseudo-time of 9.5 Gyr These curves are computed from 
a simple average of the 533 runs that produce positive matches 
with the observational data. Over the age of the Solar System, 
our results indicate that the number of Z) > 1 km Main Belt as- 
teroids striking Ceres is 4 63 1 and 1 096 for Vesta. Moreover, the 
largest Main Belt asteroids expected to have impacted Ceres and 
Vesta over the history of the Solar System have had diameters of 
71.7 and 21.1 km. 
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Fig. 3. Incremental number of Main Belt impactors onto Ceres 
(solid line) and Vesta (dashed line) as a function of impactor 
diameter obtained using a pseudo-time of 9.5 Gyr. 



3.2. Cratering on Ceres and Vesta 

Using these impactor fluxes, it is possible to obtain the number 
of craters on Ceres and Vesta as a function of crater diameter. 
To do this, we make use of Eqs. ([S) and (|4| with the values of 
the parameters Ki, v, fi, Y, and D, specified in Section 2.4, and 
the physical data shown in Table|2] Moreover, a density pi - 2.7 
g cm"^ is assumed for the impactors, which is consistent with 
the densities of several S-type asteroids (e.g. Belton et al. 1 19951 
Veverka et al. l20001 l. On the other hand, the impactor velocity Vi 
is given by 

Vi = (U^ + vUI\ (5) 

where U is the hyperbolic encounter velocity, which is equal 
to the mean impact velocity (V), and Vesc the escape velocity 
of the target. O'Brien et al. (1201 11 1 derived particular values of 
{V) equal to 4.79 km s for Ceres and 4.74 km s"' for Vesta. 
According to the physical data shown in Table|2l Ceres and Vesta 
have escape velocities of 0.51 km s"' and 0.37 km s"'. Using 
these values in Eq. (|5]), the impactor velocities used for Ceres is 
4.82 km s"' and 4.75 km s"' for Vesta. 

Figure |4] shows the total number of craters on Ceres and 
Vesta produced by Main Belt asteroids as a function of crater 
diameter obtained using a pseudo-time of 9.5 Gyr. Moreover, 
Table [3]presents a quantitative summary of our results concern- 
ing the craters formed on Ceres and Vesta. From this, the sur- 
faces of these objects should show impact structures with a wide 
range of sizes. On the one hand, our results indicate that the num- 
ber of D > 0.1 km craters on Ceres and Vesta are ~ 3.4 x 10** 
and 6.2 x 10^. On the other hand, our simulations show that 47 
craters with D > 100 km are present on Ceres and 8 on Vesta. 

Hubble Space Telescope (HST) observations of Vesta have 
revealed a singular crater with a diameter of about 460 km near 
the south pole (Thomas et al. ll997ab . Using a two-dimensional 
hydrocode, Asphaug ( 1997) studied the impact origin of Vesta 
family and concluded that the 460 km crater was formed by an 
impactor of ~ 42 km in diameter at 5.4 km s"' . Recently, Ivanov 
et al. (1201 11 1 have analyzed the formation of the south pole im- 
pact crater on Vesta from projectiles with an impact velocity of 
5.5 km s"' and sizes ranging from 40 to 96 km in diameter. Their 
results suggest the possibility that the projectile diameter may be 
larger than ~ 40 km. 
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Fig. 4. Incremental number of craters on Ceres (solid line) and 
Vesta (dashed line) produced by Main Belt asteroids as a func- 
tion of crater diameter, obtained using a pseudo-time of 9.5 Gyr. 

Table 3. Craters on Ceres and Vesta. Total number of D > 0.1, 
1,10, and 100 km craters on Ceres and Vesta produced by Main 
Belt asteroids using a pseudo-time of 9.5 Gyr. 
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Using Eqs. (|3]l from Holsapple & Housen (120071 1 and (|4| 
from McKinnon et al. (2003), the transition diameter specified 
for Vesta in Section 2.4, and the physical data shown in Table|2] 
we conclude that a D - 42 km impactor produces a transient 
crater of D ~ 217 km and a final complex crater of D ~ 307 
km on Vestifl From this, the transient crater diameter is ~ 0.7 of 
the final crater diameter. Several works such as those by Grieve 
et al. (jl98D, Melosh ( 198^, and Croft ( 1985), suggest that the 
transient crater diameter is estimated at 0.5 to 0.65 of the final 
crater diameter. However, these studies are based on lunar and 
terrestrial impact structures so we think that it would not be ap- 
propriate to compare them with our analysis relative to Vesta. 

Taking Eq. (HJl from McKinnon et al. (i200 3) into account 
for complex craters, a final crater with a diameter of ~ 460 km 
is formed by a D =42 km projectile if the transition diameter 
D, is of 2.2 km, which is a typical value for Venus and Mars. 
However, from Schenk et al. (I2004I I. we consider this value of 
the transition diameter to be too small for Vesta. Either way, the 
Dawn Mission data will be crucial for estimating the transition 
diameter for asteroid-sized objects. 

Using the most appropriate value of 50 km for the transition 
diameter D,, our study indicates that the 460 km crater observed 
on Vesta's surface had to be formed by a projectile of ~ 66.2 km 
in diameter, which is consistent with suggesting by Ivanov et al. 
(1201 11 1. Our results indicate that such an event has a probability 
of occurring of ~ 30 % over the Solar System history. We discuss 
this in the next section. 



' It is worth noting that, since the D = 460 km singular crater is com- 
parable to the Vesta diameter, the cratering laws used here may break 
down in this regime. 
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4. Discussion 

Bottke et al. (12005 al ) created a numerical model capable of track- 
ing how the Main Belt population was affected by coUisional 
evolution from the end of accretion among D < 1000 km objects 
to the present time. These authors find that the net collisional 
activity in the Main Belt over its lifetime is the equivalent of ~ 
7.5 - 9.5 Gyr of collisional activity in the current Main Belt. In 
fact, the extra comminution had to come from a collisional phase 
occurring early in Solar System history when there were many 
more D < 1000 km objects in the primordial Main Belt than in 
the current population. Finally, the extra material in this popula- 
tion would have been primarily removed by dynamical processes 
rather than collisional evolution. 

It is worth noting that there are some caveats and limitations 
to the model that should be mentioned. Bottke et al. ( I2005al) as- 
sume that most Main Belt comminution occurred when the in- 
trinsic collision probability (Pi) and the mean impact velocity 
(V) were comparable to their current values. This result is based 
on the dynamical simulations of the early Main Belt from Petit 
et al. (120011) . These authors developed a dynamical evolution 
model of the Main Belt, with Jupiter and Saturn on their cur- 
rent orbits. With this orbital configuration for the gas giant plan- 
ets. Petit et al. ( 1200 11 1 show that the depletion of asteroid Main 
Belt was very strong, removing ~ 90 % of the original popu- 
lation in less than 10 Myr after the formation of Jupiter. Based 
on this model, Bottke et al. (I2005al) assume that the dynamical 
removal phase of the Main Belt was short enough that its evolu- 
tion has been dominated by the same collisional parameters as 
found in the current Main Belt. Later, Bottke et al. (i2005b ) com- 
bined dynamical results from Petit et al. (2001) and the colli- 
sional code of Bottke et al. (2005a) to model the evolution of the 
Main Belt over the age of the Solar System. Unlike Bottke et al. 
( I2005al l. the algorithm developed by Bottke et al. (|2005b ) incor- 
porates time-varying collisional parameters such as the intrin- 
sic collision probability (Pi) and the mean impact velocity (V). 
From this model, Bottke et al. (2005b) find that the number of 
breakups occurring over 4.5 Gyr of evolution is consistent with 
the results of Bottke et al. (I2005al) using the pseudo-time approx- 
imation. Moreover, they find that the high velocity impacts do 
not appear to have a dominant effect on Main Belt comminution, 
mainly because these velocities become important only after the 
excited population has been significantly depleted. 

However, the depletion rate of the early asteroid Main Belt 
depends on the initial orbital configuration of the giant planets. 
The Nice model from Tsiganis et al. (I2005I I strongly suggests 
that the initial system of outer planets had a more compact con- 
figuration (all within ~ 15 AU of the Sun), with nearly circular 
and coplanar orbits. In the context of this model, O'Brien et al. 
( I2007I I studied the primordial excitation and clearing of the aster- 
oid Main Belt, assuming Jupiter and Saturn on initially circular 
orbits. A relevant result is that the excitation and depletion of 
the asteroid Main Belt derived from the Nice model are slower 
than those found by Petit et al. (2001) with Jupiter and Saturn 
on their current orbits. In fact, 90 % of the original population is 
removed from the Main Belt in ~ 70 Myr. Since in this scenario 
the massive population of the primitive Main Belt is removed 
more slowly, the degree of collisional evolution may have been 
greater than derived from the model of Bottke et al. (2005a). 
Therefore, the degree of collisional evolution of Ceres and Vesta 
could provide clues to the initial orbital configuration of the gas 
giant planets. 

On the other hand, the model from Bottke et al. (2005a) does 
not include the collisional effects of external impactors not be- 



longing to the Main Belt. The so-called late heavy bombardment 
(LHB) is a period of timeh during which it is believed that a large 
number of objects from the outer Solar System impacted the ter- 
restrial planets and the Main Belt asteroids about 3.9 Gyr ago. 
This event may have led to significant effects on the degree of 
collisional evolution suffered by the Main Belt. It is worth not- 
ing that the existence of the LHB is still being debated in the 
literature (Chapman et al. l2007l l. 

In particular, we estimate that Vesta had a ~ 30 % proba- 
bility of being impacted by a D = 66.2 km projectile, which is 
capable of producing the singular crater of 460 km in diameter 
observed on its surface. A slower asteroid removal rate slower 
that is consistent with the Nice model and/or the existence of the 
LHB could lead to a higher degree of collisional evolution for 
the Main Belt increasing then the probability of producing the 
crater of 460 km on Vesta. However, a ~ 30 % probability of 
occurrence does not allow us to discard the production of this 
crater in the present collisional evolution model. If the D - 460 
km crater on Vesta formed during the early massive phase of 
the Main Belt, this would be consistent with studies developed 
by Bogard & Garrison ( 120031 1. who propose that Vesta suffered 
such huge impact 4.48 Gyr ago. 

According to the comments made in this section, we believe 
that, if significant discrepancies between our results concerning 
the collisional history of Ceres and Vesta and the observational 
data obtained from the Dawn Mission were found, they should 
be linked to a higher degree of collisional evolution in the Main 
Belt. An increase in the collisional activity may be due to 

1 - a different initial configuration of the giant planets consistent 

with, for example, the Nice model, and/or 
2- the existence of the LHB. 

If such discrepancies are eventually substantiated, future works 
should focus on accurately quantifying the degree of collisional 
evolution suffered by the early massive Main Belt as well as dur- 
ing the LHB, assuming that the outer planets had initial orbits 
similar to and substantially different from (e.g. Nice model) their 
current ones. From this, the Dawn Mission would be able to give 
us interesting clues about the initial configuration of the early 
Solar System and its subsequent dynamical evolution. 

5. Conclusions 

We have presented a study aimed at analyzing the impactor flux 
and cratering on Ceres and Vesta produced by Main Belt aster- 
oids. To do this, we constructed a statistical code based on Bottke 
et al. (2005a 2005b) that includes catastrophic collisions and 
noncollisional removal processes, such as the Yarkovsky effect 
and the orbital resonances. Assuming that the dynamical deple- 
tion of the early Main Belt was very strong, and owing to that, 
most Main Belt comminution occurred when the intrinsic col- 
lision probability (Pi) and the mean impact velocity (V) were 
comparable to their current values. Our main results are the fol- 
lowing 

- The number of D > 1 km Main Belt asteroids striking Ceres 
and Vesta over the Solar System history is approximately 
4600 and 1 100. 

- The largest Main Belt asteroids expected to have impacted 
Ceres and Vesta had diameters of 71.7 and 21.1 km. 

- The surfaces of Ceres and Vesta should show a wide variety 
of craters with a wide range of sizes. On the one hand, the 
number of D > 0.1 km craters on Ceres and Vesta are ~ 3.4 
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X 10^ and 6.2 x 10^. On the other hand, the surfaces of Ceres 
and Vesta present 47 and 8 craters with D > 100 km. 

- Using an appropriate value of 50 km for the transition diam- 
eter from simple to complex craters together with the expres- 
sions derived by Holsapple & Housen 12007 ) and McKinnon 
et al. (I2003I I. the D = 460 km crater observed on Vesta had 
to be formed by a D ~ 66.2 km projectile. 

- Such an event has an occurrence probability ~ 30 % over 
the Solar System history. Then, we cannot discard the pro- 
duction of this crater in the present collisional evolution 
model. An asteroid removal rate that is slower is consistent 
with the Nice model and/or the existence of the LHB which 
could lead to a higher degree of collisional evolution of the 
Main Belt thereby increasing the probability of producing 
the crater of 460 km on Vesta. 

- We suggest that, although the impact that formed the D - 
460 km crater on Vesta cannot be discarded in the present 
collisional evolution model, it is more likely to have occurred 
in the early massive phase of Main Belt while it was be- 
ing excited and depleted or during the LHB. This conclusion 
may become consistent with studies from Bogard & Garrison 
(I2003I I who propose that such an impact occurred 4.48 Gyr 
ago. 

- We believe that if significant discrepancies between our re- 
sults about the cratering on Ceres and Vesta and data ob- 
tained from the Dawn Mission were found, they should be 
linked to a higher degree of collisional evolution during the 
early massive phase of the Main Belt and/or the existence of 
the LHB. From this, the Dawn Mission could play an im- 
portant role by providing information on the formation and 
evolution of the early Solar System. 
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